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ABSTRACT 

Context. The distance to an isolated dark globule is often unknown and yet crucial for understanding its properties, in particular its 
mass. A new approach to this problem is discussed 

Aims. The purpose of the present paper is to investigate how well the distances of more or less reddened field stars can be determined 
by using multi-colour imaging. 

Methods. We observed a test globule, B 335 in U, B, g, r, and I, and together with the 2MASS survey, this data set gives a well-defined 
spectral energy distribution (SED) of a large number of stars. The SED of each star depends on the interstellar extinction, the distance 
to the star, and its intrinsic SED. As we had good reasons to suspect that the wavelength dependence of the extinction (the reddening) 
changes from the outskirts of the globule to the central parts, we did not assume any specific reddening law. Instead, we use a scheme 
that allows independent determination of the extinction in each line of sight as determined by groups of adjacent stars. The method is 
based on the use of stellar atmospheric models to represent the intrinsic SEDs of the stars. Formally, it is then possible to determine 
the spectral class of each star and thereby its distance. For some of the stars we have optical spectra, allowing us to compare the 
photometric classification to the spectrometric. 

Results. As expected, the main problem is that there are few stars found within each distance bin for the small field size defining a 
typical dark globule. However, the characterisation of the extinction and photometric classification give consistent results and we can 
identify one star at the front side of the globule. It has a photometric distance of 90 pc. The closest star behind the B 335 globule has 
a distance of only a 120 pc and we therefore determine the distance to B335 as 90-120 pc. Our deep U image shows a relatively 
bright south-western rim of the globule, and we investigate whether it might be due to a local enhancement of the radiation field. A 
candidate source, located 1.5 arcminutes outside our field, would be the field star, HD 184982. This star has an entry in the Hipparcos 
Catalogue and its distance is 140-200 pc. However, we come to the conclusion that the bright SW rim is more likely due to the wing 
of the point-spread-function (PSF) of this star. 

Key words. ISM: clouds - ISM: individual objects: B 335 - ISM: dust, extinction 



1. Introduction 

It is in principle easy to determine the distance to a dark cloud: 
it is just a matter of finding a star just in front of the cloud and 
another closely behind it. Or, even better, a star close enough to 
the cloud to give rise to a reflection nebula. The closest cloud 
complexes, like that in Taurus/Aurigae, cover large regions on 
, the sky and thanks to the Hipparcos mission ( 1997 J there are a 
■ number of stars with known distances within the field. The re- 
' maining task is that of placing the stars in front of or behind 
, the cloud. The traditional way is to look for the onset of inter- 
stellar reddening (see Lombardi et al. 120081 for a recent exam- 
ple). As an attractive alternative, the onset of polarisation can be 
determined (Alves & Franco l2007l l. For small isolated clouds, 
however, there are too few (if any) Hipparcos stars within the 
field, and the distances to the stars must be determined by spec- 
troscopic/photometric means. The reason why it is important to 
determine the distances to such clouds, including dark globules, 
is their role of hosting isolated star formation. Thus, important 
properties as the luminosity of the forming star and the mass and 
the radius of the cloud cannot be determined without knowing 
the distance. 



* Based on observations collected at the European Southern 
Observatory, Chile (ESO programme 077.C-0524(A)) 



Maheswar & Bhatt (I2006I I have estimated the distances to 
nine dark globules using V, R and I photometry in combina- 
tion with the J, H and K photometry from the 2MASS catalogue 
(Skrutskie et al. 2006). By assuming a standard interstellar ex- 
tinction law (valid for diffuse clouds) and adjusting visual ex- 
tinction. Ay, until the SED fits that of a main sequence (MS) 
star, these authors determine the extinction and the distance for 
each star. There are two uncertain assumptions in their approach, 
which may have strong impact on the estimates of the distances. 
The first is the assumption of a standard extinction law, while 
it is known that the Ry value (/?i/=Av/E(B-V)) changes from 
typically Rv=3.1 in the diffuse ISM to Ry ~ 5 in molecular 
clouds (see e g Ossenkopf & Henning 1994 Strafella et al l2001l l. 
The second is the assumption that all the stars are MS stars. In 
the present paper we discuss these problems and present a more 
elaborate method which is based on imaging in eight filters. 

As a test case we use the well studied dark globule B 335. It 
hosts a protostar heavily obscured by a flattened core structure, 
including an accretion disk seen edge-on (Galfalk & Olofsson 
2007) . In most of the numerous investigations of this source, a 
distance of 250 pc as given by Tomita et al. ( 1979^ is assumed. 
Stutz et al. (2008) have recently revised the distance estimate by 
using archive photometry of stars within a radius of 20'. Their 
revised distance is 150pc (60-200pc). The diameter of B335 
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" co-added image provided by M. Galfalk 



as seen in the optical is just » 4' . CO observations (Frerking et 
al.[l987 ), however, show that it extends far to the NE, including 
stars within such a large radius as 20'. 

2. Observations and data reduction 

B335 (RA(2000) = 19 37 0.8, DEC(2000) = 7 34 11 ]) 
and a reference field (centered at 19 34 10., 07 37 05) were 
observed using the ESO/NTT at La Silla during four nights 
2006-06-27-30. The reference field was chosen to have the 
same galactic latitude as B 335. The observational details are 
given in Table [1] The EMMI images have been complemented 
by a Kj image taken with the 2.5 m Nordic Telescope on La 
Palma (total exposing 48 min; cf. Galfalk et al. 2007 ). 
Furthermore, spectra from ten of the stars in the field were 
measured with NTT-EMMI in RILD-mode with a resolution of 
~ 600 in the wavelength range 550 - 1000 nm (exposure times 
between 300 and 1000 s, sUt 1.5" and 5 "). 
All object exposures were preceded and followed by standard 
star (Landholt equatorial standards SAl 11-2093) observations 
in respective filter. The basic reductions (bias subtraction, dark 
coiTection, flat-fielding, cosmic ray reduction) were carried 
out using standard IRAF routines. Then the IRAF astro me trie 
programs were used for registering and co-adding the images 
using Sky view (http://skyview.gsfc.nasa.gov) facilities for 
astrometric data on stars common to all images of an object. 
The summation of the frames were weighted by the inverse of 
the errors of the stars in the middle magnitude range. The star 
finding program 'sextractor' were then used to tabulate stars 
and positions. Photometry on all co-added images were then 
performed using the DAOPHOT photometry package. 
In addition, our observations were combined with data from the 
2MASS Point Source Catalog (PSC) (Skrutskie et al. 2006) thus 
extending the stellar data with the J-H-K measurements (where 
these exist) in the B 335 and the reference field. Furthermore, 
object field data in wavelengths around 3.6, 4.5, 5.8, 8.0;U from 
the Spitzer data archive were used. The Spitzer archived images 
were analyzed with the Mopex-software package (see Makovoz 
& Kh an l2005l Makovoz & Marleau 1^051 and Masci et al 
12005 1) to give the stellar magnitudes. Finally, tables for B 335 
and the reference field containing object positions, magnitudes 
and magnitude errors in all filters were compiled. The field of 
view differs between the red and the blue arm of EMMI and 
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Fig. 1. Unreddened U-I and I-K colours for giants (full lines) and 
main sequence stars (dashed) as a function of effective tempera- 
ture based on the model atmospheres of Hauschildt et al.(1999). 
K and M giants can be separated from main sequence stars with 
the same temperature by measuring the U-I colour index. The I 
- K colour index is sensitive to gravity for late M stars. 



the 4'field of NOTCam was not exactly centered at the same 
position as the EMMI field. This results in the observation field 
which is slightly less than 6 arc-minutes. 



3. Results and interpretation 

3.1. Setting the scene 

Our observations basically gave the SEDs for a large number of 
stars, more or less hidden by the globule. Our purpose was to 
identify stars immediately behind and in front of the globule - 
if there are any. To do so, we needed the intrinsic colours of the 
different stars and their corresponding absolute magnitudes. The 
single most important parameter that determines the colour of a 
star is its temperature, but for a given temperature the luminosity 
is very different for a dwarf as compared to a giant or supergiant 
star supergiants are so rare that the chance of finding one within 
such a small field, a few degrees away from the Galactic plane, 
is small. As will be discussed below, giant stars cannot be ex- 
cluded on the same grounds, and we had to look for means to 
separate giants from dwarfs. To avoid the difficulties to inter- 
polate between different photometric systems, we chose to use 
synthetic colours to represent the intrinsic colours of the stars. 
We chose the Next Generation model atmospheres provided by 
Hauschildt et al. 119991 In order to keep the parameters within 
reasonable limits, we restricted the range of surface gravity to 
one value typical for MS stars (logg = 4.5) and one typical for 
giants (logg - 0). We also restricted the metal abundance to so- 
lar (typical for Population I, which dominates in most directions 
along the Galactic plane). As one may suspect, it is hard to sepa- 
rate dwarfs from giants in most combinations of colour indices. 
However, there is a tendency for the cool (K and M) giants to be 
redder than the dwarfs for the same T^fj and this effect is pro- 
nounced in the U-I colour index, but also seen in eg. the I-K 
index for late M stars (see Fig [T]l. By combining these two in- 
dices in a colour-colour plot, we should be able to separate K 
and M giants from dwarfs (Fig|2]i. 
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Fig. 2. The (U-I) colour versus (I-K) colour plot has the potential 
to separate K and M giants from main sequence stars - provided 
the interstellar reddening vector E(U-I)/E(I-K) is less than the 
slope of the straight part of the curves 2.7). 
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Fig. 3. The accumulated number of stars within 25 square ar- 
cminutes as a function of distance as predicted by the galaxy 
model of Wainscoat & Cohen (il992i) for the direction of B 335. 

In a volume limited sample, dwarfs are much more frequent 
than giants which means that it is very likely that the closest 
stars are dwarf stars. On the other hand, for a given apparent 
flux limit, the horizon is much further away (typically an order 
of magnitude) for giants than for dwarfs, and giants will be over- 
represented by as much as three orders of magnitude. In order to 
quantitatively estimate this effect, we used the stellar distribu- 
tion model of the Galaxy as described by Wainscoat and Cohen 
( fT992] l for the galactic coordinates of B 335 (/=44.9, b^-6.6). In 
Fig|3]we show the expected number of stars within our field of 
view as a function of distance. As expected, the MS stars domi- 
nate. On the other hand, the giants dominate among the brightest 
stars (FigHJi. 
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Fig. 4. The K luminosity function for a 25 square arcminutes 
field as predicted by the Galaxy model of Wainscoat & Cohen 
( (19921 ) for the direction of B 335. 

3.2. How fo identify tiie closest stars? 

In Fig|5]we display an infrared HR diagram with I-K as colour 
(temperature) axis and the apparent K magnitude as the bright- 
ness axis. We also indicate the MS locus for three different 
distances. The brightest star in the field is a known variable 
star (V1142 Aql, equivalent to S 9453 in the catalogue of 
Hoffmeister 1 1966) and it is probably an AGB star. Because of 
its extreme position in the diagram (I-K = 13.6), we left it out- 
side the plotting range. The positions of the observed stars in the 
diagram is determined by the luminosity and the colour as well 
as the extinction and the distance for each star. We have also in- 
dicated the reddening vector for the K-extinction in the diagram. 
The slope of the reddening vector is less than that of the locus 
for the main-sequence stars in the diagram. Thus correcting for 
the extinction will result in still longer distance estimates. In ad- 
dition, as discussed above, most of the stars seemingly within 
150 pc may well be giants far behind the cloud. 

From Fig|6]we could conclude that all the 'nearby' stars 
appear to be giants, except one which is a heavily obscured 
early-type star. All of these are therefore remote, except pos- 
sibly No. 2. With the following reasoning we find that this star 
also is remote. The reddening line through star No. 2 crosses the 
giant branch at two positions with T^ff ~ 3400 K and a; 4200 K 
(which corresponds to M5III and K4III respectively). If it is a 
M5III star, the reddening is E(I-K) - 2.0 which correspond to 
Ak ~0.5. The apparent magnitude of star No. 2 is K=7.51 and 
as the absolute magnitude of a M5III is Mk ~ -5, the distance 
to this star would be x3 kpc. 

We next turned to the 23 stars seemingly at a distance be- 
tween 75 and 150pc and display them in a U-I versus I-K dia- 
gram (Fig|7]i. Except for the four stars closest to the lower red- 
dening line (which probably are early-type MS stars), it is not 
certain whether to designate the stars as dwarfs or giants be- 
cause the reddening vectors intersect both the locus of the giants 
and that of the dwarfs. Although the general impression is that 
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the stars are heavily obscured (and therefore behind the cloud), 
we could not exclude that a star close to the cool part of the MS 
locus actually is a foreground star. If this were the case, its dis- 
tance would represent a lower limit of the distance to the cloud. 
To further constrain the properties of these stars we obviously 
needed more information and as we, so far, only had considered 
one colour-colour diagram out of several possible it is natural to 
look for some other combination. It is well known that the con- 
tinuum opacity minimum at 1 .6 i-im is pressure dependent, which 
is reflected in the J-H versus H-K diagram (Fig|8]l. We marked 
the position of the four possible foreground stars in this diagram, 
and indeed, one of them (No. 12) can either be a late M dwarf or 
(if behind the cloud) a star with Te/f ~ 6000 K. However, such 
a high temperature can be excluded from the position of this star 
in the U-I versus I-K diagram, and this star is therefore very 
likely a cool foreground star. 

3.3. A more efficient approacli 

In principle it would be possible to further constrain the proper- 
ties of the stars by inspection of other colour combinations, but 
it is not a very practical way to proceed. In addition, we have so 
far assumed that the reddening law is independent of the light 
path through the cloud, be it close to the center or in the out- 
skirts. This is not likely to be the real case but, as we will show 
in a forthcoming paper, the scheme introduced by Cardelli et 
al.(Il989) where the extinction is defined by a single parameter 
Rv, adequately represents the extinction. 

A more general and eflicient approach would be to start by 
searching for stars that are likely to be MS stars. With this step 
we sought the best ;^'-square fit of the SED for each star among 
the 61 model SEDs, optimising for the Rv parameter and the ex- 
tinction (Ay). This optimisation is fast and simple, but one may 
wonder whether it is reasonable to let the Ry parameter be a free 
parameter for each star. Actually, when correlating the resulting 
Rv parameters to the positions of the stars, we saw a tendency in 
the sense that the /?v-values are higher in the eastern part than in 
the western. This will be discussed further in a forthcoming pa- 
per, but more significant for the present analysis is that we noted 
a high scatter regarding the /?v-values for adjacent stars. This is 
not likely to reflect the true conditions, and therefore we intro- 
duced a second step. In this step we selected probable nearby 
MS stars found in the first step. For each of them we repeat the 
optimisation, but this time together with a few adjacent stars re- 
questing that this grouping experiences a cloud extinction char- 
acterised by the same Ry-value. We allowed, however, for an 
individual total extinction, A[/. Finally, in order to get a measure 
on the stability of the solutions we repeated the calculations a 
few times for diff'erent combinations of adjacent stars. The re- 
sults are shown in Fig|9] We first noted that star No. 12 also in 
this analysis is found to be a foreground star, and its photometric 
distance is 90 pc. The two closest stars behind the cloud is No. 9 
(120 pc) and No. 13 (140 pc). The distance to B 335 is therefore 
90-120pc, neglecting possible systematic errors. 

3.4. Can we trust tiie pliotometric ciassification? 

We classified the collected spectra by comparison to templates 
from the ESO/UVES survey (12005]) . The temperature, T^//, 
is then taken from the compilation in Allen's Astrophysical 
Quantities (2000). The photometric classification gives T^ff and 
the preferred status giant-or-dwarf, and for convenience we used 
the same calibration to indicate the spectral class. The result is 
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Fig. 5. In this infrared HR diagram we have plotted the stars in 
the B 335 field. The curves represent the main sequence stars at 
75, 150 and 300 pc. The reddening vector indicates the effect of 
increased extinction on a star position. The seemingly closest 
stars are identified by larger symbols and numbers. 

shown in Table|2] The agreement is in most cases good. 
It is, however, worrying that the only star in the list that may be 
of interest for the distance estimate. No 13, differs considerably. 
The TiO and CaH bands are pronounced so it is definitely a M 
star and the Nal lines indicate that it is a dwarf star. However, 
it cannot be completely excluded that it may be a giant and that 
the strength of the sodium doublet is due to interstellar absorp- 
tion components. In either case, this star is not relevant for our 
distance determination, as it not the closest star behind the cloud. 



3.5. The briglit rim of B 335 

Dark globules often show bright rims due to scattered light orig- 
inating from the surrounding illumination by the interstellar ra- 
diation field. One may wonder why the core of the globule looks 
dark as it is also illuminated by the interstellar radiation field. 
The reason is the angular dependence of the scattering, being 
more efficient in the forward direction. To illustrate this point 
we calculated the phase diagram for the grain size distribution 
model (MRN) given by Mathis et al (1977). 

In Fig [To] we show the scattering efficiency as a function of 
angle, and we include both the wavelength for the U and the B 
band. The forward scattering is x factor of 4 larger than the back 
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Fig. 6. In this colour-colour diagram we have only included 
stars with good photometry (cr < 0.05 magnitudes). The seem- 
ingly closest stars are marked by larger symbols (cf. FigjSj. The 
dashed curve denotes the the expected locus for giants and the 
full drawn curve the locus for MS stars (No. 20 being one of 
them). The lowest temperature is 3000 K, and each 200 K is in- 
dicated. The reddening vector is well defined by the locus of 
early type MS stars. All the seemingly closest stars, except No. 
20, are most probably giants. 
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Fig. 7. In this colour-colour diagram including the stars within 
the 75 and 150 pc curves in Fig |5] there are no obvious giants. 
On the other hand, tracing back along the reddening vector gives 
intersections with both the giant and MS curves, except possi- 
bly for the five stars closest to the lower reddening line, which 
appears to be early type MS stars. The four encircled stars are 
potential candidates for foreground M dwarfs. 



scattering, explaining why the rim of a dark globule in general 
is bright. It must also be realized that it is virtually impossible 
to determine the sky emission accurately enough to measure the 
light level at the center of a dark globule. In the case of B 335, we 
do observe a bright rim and we note a pronounced enhancement 
in the south-west indicating that there is a local enhancement due 
to a nearby star. We first considered the case that this star illu- 
minates the globule. We next investigated the possibility that the 
light enhancement is due to the point spread function (PSF) of a 



Fig. 8. In this J - H colour vs H - K colour diagram the four stars 
identified in Fig|7]as possible foreground stars are marked with 
larger symbols. We argue that one of them. No. 12, is indeed a 
foreground M dwarf. 
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Fig. 9. This diagram shows the distances and the extinctions 
to the closest stars in the B 335 field. Only one star is found 
in front of the cloud. The error bars only reflect the ;t'-square 
fitting errors, and the indicated distance to B 335 is 90-120pc. 
The dashed line schematically represents the extinction due to 
the diff^use ISM. 



bright star located outside our field. The most probable star (in 
both cases) is HD 184982 which is located 1.5 arcminutes to the 
south-west of the image. This star has an entry in the Hipparcos 
Catalogue and its distance is 140-200 pc. 

As the center of the globule is dark, this star - if physically 
associated - must be located behind the cloud, but how far be- 
hind? We first note that the surface brightness close to the south- 
western edge of the image, U = 24 magnitudes arcsec"^, corre- 
sponds to 0.4 yuJy arccsec which is comparable to the inter- 
stellar radiation field, l.OyuJy arcsec"^. We must, however, keep 
in mind that the rim of the globule only is illuminated from one 
side and that the albedo of the grains is less than unity. Therefore, 
additional illumination is required to explain the bright south- 
western rim of the globule. The A2V star has an absolute mag- 
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Fig. 10. The phase diagram for the MRN grain size distribution 
model. It is normahsed for the B filter at a scattering angle of 
zero degrees. The diagram shows that the forward scattering is 
much more efficient than the back scattering and also, that this 
is more pronounced in the U than in B band. 



nitude M(U)^1.4 which means that it gives the same (one-side) 
illumination as the interstellar radiation field if it were located 
0.4 pc from the cloud. It is likely that the star is located a bit 
closer to the cloud in order to give rise to the observed illumina- 
tion, but for our purpose it suffices to conclude that the distance 
to the B 335 globule is essentially the same as that to HD 184982 
(provided that it really is physically associated). We now turned 
to the possibility that we just measured the wing of the point 
spread function (PSF) of this star. By co-adding 23 images of 
the calibration star, we are able to define the PSF to a radius 
of 40 arcseconds. A power-law extrapolation showed that the 
observed rim has a surface brightness of a factor 2 higher than 
expected from the PSF (Fig[T2li. The index of the power law ex- 
trapolation, -2.6, makes sense, but in view of how far (from a 
radius of 40 to 100 arcseconds) we had extrapolated the mea- 
sured PFS, we could not claim that the difference between the 
measured surface brightness at the SW rim of B 335 is signifi- 
cantly higher than the PSF wing of HD 184982. Therefore this 
explanation is (unfortunately) more likely than the physical as- 
sociation. If the photometric distance to the globule of 90-120 pc 
is accepted the star HD 184982 (with its Stromgren photometric 
distance of 127pc) is still one of the stars close to the cloud. 

4. Discussion 

4. 7 . How fo apply the method to a large sample of globules? 

If the Galactic stellar distribution model by Wainscoat & Cohen 
( 11992 1) applies, then we were very lucky to find four stars in our 
field being closer than 200 pc. The model may underestimate the 
low-mass population of stars, and if so, our field of view (a;36 
square arcminutes) may suffice. Still, in planning for a survey 
of dark globules aiming at distance determinations, it would be 
better to observe larger fields. However, in most cases, the glob- 
ules as seen in the optical are small and it is often difficult to 
define how far their halos extend. Even though the UV imaging 
is more efficient in revealing the extent of the globule (includ- 
ing its bright rim), a separate CO survey designed to define the 
sizes would be very useful. In our pilot study we spent a very 




Fig. 11. The U image of B 335 shows a rim of scattered light as 
is often the case for dark globules. The pronounced enhancement 
in the S W corner is proposed to be due to a bright star 1 '.5 outside 
the field. The lowest contour level corresponds to U = 27 arcsec^ 
and the step size is 0.5 magn 



long time on the UV imaging. This was justified by the expecta- 
tion that the closest stars would likely be M dwarfs and therefore 
faint in the UV and, in addition, the extinction of the globule is 
higher at short wavelengths. The faintest of the nine closest stars. 
No 35, has a UV brightness of U = 20.00 (cr=0.014), and as a 
photometric accuracy of cr=0.04 is estimated to be good enough, 
1 h exposure time should have sufficed (instead of the 8 h spent). 
Also for the other photometric bands in the optical, the exposure 
times could have been much shorter. The faintest star in the Kj 
band among these nine closest star. No 25, is bright, K, = 1 1.75 
(cr=0.002), and consequently the 2MASS survey provides the 
required accuracy. This all means that typically 2-3 hours per 
globule should suffice for a distance survey using a wide-field 
camera on a medium sized telescope. 

4.2. is the method reliable? 

It is clear from this test case that a larger field than our 36 square 
arc-minutes should be used in order to include more stars at 
the relevant distances. The photometric classification is in most 
cases adequate, but we noted exceptions. One reason for ex- 
ceptions is probably significantly lower metal abundance. It is 
known (Buser & Kurucz, 119921 ) that colour indices including 
the U band do depend on the metal abundance. In order to quan- 
tify this effect we have used a set of model atmospheres with 
[Fe/H] =-1, and the difference compared to solar abundance in 
the U-I colour index is 0.3 magnitudes for late M stars. This is 
far more than the photometric accuracy. It is likely that the mis- 
match for star No 13, discussed above, is caused by a low metal 
abundance. Again, a larger field would include more stars at a 
relevant distance and make it easier to exclude stars with very 
different metal abundance compared to the majority of the stars. 
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and as the cool MS stars are the most abundant, any method 
should be able to sort out these stars and reliably separate 
them from cool giant stars. 

2. It turns out that multi-colour imaging in combination with 
synthetic colours can provide reliable spectral classification 
and extinction, and we find that the U band is particularly 
important, both for the spectral classification and for tracing 
the extent of the globule. 

3. We have tested a new SED-fitting method which allows for 
different Rv-values across the field. For the test case, B 335, 
we get a distance of 90-120 pc. 

4. Our deep U image of B 335 shows a bright halo due to scat- 
tered star light, and the south-western rim is much brighter. 
We find that a possible explanation is that the globule is illu- 
minated by a star, HD 184982, located 1.5 arcminutes to the 
south-west of our field and a few tenths of a parcec behind 
the cloud. The distance to this star is 140-200pc according 
to the Hipparcos Catalogue. However, the surface brightness 
can also be explained as due to the PSF wing of the same star, 
and we judge that this is in fact a more hkely explanation. 




30 I I 

0.0 0.5 1.0 1.5 2.0 2.5 

log r (arcsec) 

Fig. 12. The point-spread-function of HD 184982. The extrapo- 
lated PSF wing falls sUghtly below the measured surface bright- 
ness at the south-western rim of B 335 (crosses). 
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In any case, follow-up spectroscopic observations of the stars 
identified by the photometry as being located closely in front 
and behind the cloud would be useful for verifying (refining) the 
spectral classification. 

One obvious source of uncertainty is non-resolved binaries. 
Thus, if the stars have approximately the same spectral class, the 
inferred distance would be too small by a factor V2. Otherwise 
the error will of course be less, and there is also a chance that our 
SED-fitting procedure may flag such cases (simply by a worse 
fit), but we have not looked closer at this aspect in our pilot study. 
But, again, this compUcation is less of a problem when a larger 
field is applicable improving the statistics. 

4.3. Is the method computationally practical? 

The first step of our method, i.e. finding the best SED fit with 
the distance and the Ry-value as parameters, is fast and simple. 
The reason why we found a need for a second step is the result- 
ing large scatter of the Ry-value between adjacent stars, which 
is not likely to be real. To cope with this problem, we added the 
second step where we carried out the SED fitting for adjacent 
stars requesting that they share the same Ry-value. The compu- 
tational time for this step is proportional to the Nth power, where 
N is the number of stars in each group. In order to get a measure 
on the stability of a fit, we also ran the optimisation for different 
combinations of adjacent stars. In practice we found that N - 3 
is sufficient and that the computational time can be kept low. 
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5. Conclusions 



1. The fundamental limitation when it comes to determining 
the distance to a small dark cloud is the space density of stars 



